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Summary. After a brief introduction to the Painleve property for ordinary differ- 
ential equations, we present a concise review of the various methods of singularity 
analysis which are commonly referred to as Painleve tests. The tests are applied 
to several different examples, and the connection between singularity structure and 
integrability for ordinary and partial differential equation is discussed. 



1 Introduction 

The connection between the integrability of differential equations and the 
singularity structure of their solutions was first discovered in the pioneering 
work of Kowalevski |65j , who considered the equations for the motion under 
gravity of a rigid body about a fixed point, namely 

d£ — ixLj + cxg , 
di 

dg = gxtj ; 1 = 1^- 

In the above, £ and w are respectively the angular momentum and angular 
velocity of the body, g is the gravity vector with respect to a moving frame, 
and the centre of mass vector c and inertia tensor I are both constant. The 
remarkable insight of Kowalevski was that the system of equations ([T]) could 
be solved explicitly whenever the dependent variables £ and g are meromorphic 
functions of time t extended to the complex plane, i G C. By requiring that 
the solutions should admit Laurent expansions around singular points, she 
found constraints on the constants c and I = diag(/i, /2, /a) (diagonalized in 
a suitable frame) . Her method isolated the two solvable cases previously known 
to Euler (c = 0) and Lagrange (/i = I2 with c defining the axis of symmetry) , 
as well as a new case having principal moments of inertia Ii = I2 = 2/3 and 
c perpendicular to the axis of symmetry. The latter case is now known as 
the Kowalevski top, and Kowalevski was further able to integrate it explicitly 
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in terms of theta- functions associated with a hyperelliptic curve of genus 2, 
thereby proving directly that the solutions are meromorphic functions of t. A 
modern discussion can be found in [5] or 'M], for instance. 

An important feature of the equations ([1]) from the point of view of sin- 
gularity analysis is that they are nonlinear. For a linear differential equation 

— + Q»-i(^) ^^n-i + • • ■ + "i*-^-*^ + "■oizjy = 0, 



of arbitrary order n it is well known |43[ 155] that the general solution can have 
only fixed singularities at the points in the complex z-plane where the coef- 
ficient functions aj{z) are singular. However, for nonlinear differential equa- 
tions, as well as the fixed singularities which are determined by the equation 
itself, the solutions can have movable singularities which vary with the initial 
conditions. For example, the first order nonlinear differential equation 



dy 
dz 



has the general solution 



y = , zq arbitrary, 

with a movable simple pole at z = zo- If the initial data y = yo is specified at 
the point z — 0, then the position of the simple pole varies according to 

1 

zo = . 

yo 

The classification of ordinary differential equations (ODEs) in terms of 
their singularity structure was initiated in the work of Painleve [T^- The 
main property that Painleve sought for ODEs was that their solutions should 
be single- valued around movable singular points. Nowadays this property is 
usually formulated thus: 

Definition 1.1. The Painleve property for ODEs: An ODE has the 

Painleve property if all movable singularities of all solutions are poles. 

Painleve proved that for first order ODEs of the general form 

^ Q{y,zy 

where V and Q are polynomial functions of y and analytic functions of z (and 
the prime ' denotes d/dz), then the only movable singularities that can arise 
are poles and algebraic branch points. The latter are excluded by Definition 
1.1, and he further showed that the most general first order equation with the 
Painleve property is the Riccati equation 
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V' = 02(2)2/^ + ai{z)y + 00(2;), 

where the coefBcients ap, ai, 02 are analytic functions of z. 

For second order ODEs, life is more complicated because movable essential 
singularities can arise (see e.g. chapter 3 in [6] for an example). Painleve 
initiated the classification of second order ODEs of the form 

y" = T{y',y,z). (2) 

with !F being a rational function of y' and y, and analytic in z. Painleve and 
his contemporaries succeeded in classifying all ODEs of the type ([2|) which 
fulfill the requirements of Definition 1.1. The complete result is in the form 
of a list of approximately fifty representative equations, unique up to Mobius 
transformations, which are detailed in chapter 14 of Ince's book [55]. It was 
found that (after suitable changes of variables) all of these ODEs have general 
solutions in terms of classical special functions (defined by linear equations) 
or elliptic functions, except for six special equations which are now known as 
Painleve I- VI (or just PI- VI). 

As an example, consider the second order ODE 

y" = 62/^-^52. (3) 

This can be immediately integrated once, because the equation is autonomous 
(the right hand side is independent of z), which yields 

{y'f = 4y3 _ g^y - (4) 

with (73 being an integration constant. The general solution of the first order 
ODE ^ is given by the Weierstrass elliptic function, 

y ^ p{z ~ za;g2,g3) (5) 

with the constants g2, gs being the invariants. The solution ([5]) has in- 
finitely many movable double poles, at z = zq and at all congruent points 
z = zq -f 2muji + 2nuj2 G C for (m, n) G 1? on the period lattice defined by 
the half-periods uji,uj2. (For an introduction to Weierstrass elliptic functions 
see chapter 20 in 102 or chapter VI in ^91j.) We shall return to the equation 
^ in the next section. 

The first of the Painleve equations is PI, which is a non-autonomous version 
of ^ given by 

y" = + z. (6) 

Its general solution is a meromorphic function of z, and the solution of PI (or 
sometimes the equation itself) may be referred to as a Painleve transcendent, 
since it essentially defines a new transcendental function. The other equations 
PII-PVI also contain parameters; for example the second Painleve equation 
(PII) is 
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(7) 



where a is a constant parameter. The general sohition of each of the Painleve 
equations cannot be expressed in terms of elhptic functions or other classical 
special functions [57 , although for special parameter values they can be solved 
in this way; e.g. when a is an integer, equation ([7|) has particular solutions 
given by rational functions of 2, and it has solutions in terms of Airy functions 
for half-integer values of a. 

An important early result was the connection of PVI with the isomon- 
odromic deformation of an associated linear system ,31|. After the work of 
Painleve and his colleagues around the turn of the last century, the Painleve 
equations were probably only of interest to experts on differential equations. 
However, in the latter half of the 20th century the Painleve transcendents en- 
joyed something of a renaissance when it was discovered that they gave exact 
formulae for correlation functions in solvable models of statistical mechanics 
|103j , quantum field theory and random matrix models [27l [59] , and also 
arose as symmetry reductions of partial differential equations (PDEs) inte- 
grable by the inverse scattering transform (see [3] and section 3 below). The 
link with integrable PDEs and linear Lax pairs established the exact solution 
of the Painleve equations by the isomonodromy method . More recently a 
weakened version of the Painleve property has been used to find exact metrics 
for relativistic fluids [W. With this wide variety of physical applications, the 
Painleve transcendents have acquired the status of nonlinear special functions 
(see the review and references in chapter 7 of [5]). 

The continuation of Painleve's classification programme to higher order 
equations becomes increasingly difficult as the order increases. Even at third 
order a new phenomenon can be encountered, in the form of a movable nat- 
ural barrier or boundary beyond which the solution cannot be analytically 
continued; this occurs in Chazy's equation 



A variety of results for third or higher order equations have been obtained 
by Chazy [TJ], Gambler, Bureau, and most recently by Cosgrove; see [23] 
and references therein. Chazy's equation ([8]) has some higher order relatives 
known as Darboux-Halphen systems, which have a very complicated singular- 
ity structure, and occur as reductions of the integrable self-dual Yang-Mills 
equations (see the contribution of Ablowitz et al in ^22]). 

It should be clear from the above that the Painleve property has a very 
deep connection with the concept of integrability. This connection is by no 
means straighforward, and continues to be the subject of active research 122] . 
In the rest of this brief review article, we will introduce the basic techniques 
for testing the singularity structure of differential equations (both ODEs and 
PDEs), which are often referred to collectively as Painleve analysis. The ba- 
sic method for testing ODEs by expansions in Laurent series is treated in 
section 2. This method should probably be referred to as the Kowalevski- 
Painleve test to honour both pioneers of the subject, but most commonly 



y"' = 2yy"-3(y')'. 



(8) 
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only Painleve is mentioned in this context. Section 3 describes the conjec- 
ture of Ablowitz, Ramani and Segur [1] on the connection between integrable 
PDEs and Painleve-type ODEs, and how this can be used as an integrability 
test for PDEs. In the fourth section we explain how the preceding analysis can 
be bypassed by a direct consideration of the singularity structure of a PDE, 
by using the method of Weiss, Tabor and Carnevale [lOOj . This is followed in 
section 5 by associated truncation techniques related to Backlund transforma- 
tions, Lax pairs and Hirota bilinear equations for integrable systems of PDEs. 
In section 6 we highlight the limitations of the Painleve property as a criterion 
for integrability, in the context of integrable systems with movable algebraic 
branching and the weak Painleve property [86| . In the final section we give 
our outlook on methods of singularity analysis for differential equations, and 
mention how some of these methods and concepts have been extended to the 
discrete context of maps or difference equations. 

2 Painleve analysis for ODEs 

Consider an ODE for a dependent variable y{z), which may be a single scalar, 
or a vector quantity. If the ODE has the Painleve property then its solutions 
must have local Laurent expansions around movable singularities at z — zq, 
where zq is arbitrary. However, if branching occurs then this can be detected 
by local singularity analysis. The basic Painleve test for ODEs consists of the 
following steps: 

• Step 1: Identify all possible dominant balances, i.e. all singularities of form 
y ^ co{z- zo)^- 

• Step 2: If all exponents /i are integers, find the resonances where arbitrary 
constants can appear. 

• Step 3: If all resonances are integers, check the resonance conditions in 
each Laurent expansion. 

• Conclusion: If no obstruction is found in steps 1-3 for every dominant 
balance then the Painleve test is satisfied. 

Note that the exponents /i and leading coefficients cq must have as many 
components as the vector y, and if the ODE is polynomial then at least one 
of the exponents must be a negative integer for a leading order pole-type 
singularity. Rather than give formal definitions of the terms introduced in 
steps 1-3 above (which can be found in [T3] and elsewhere in the references), 
we would like to illustrate them with a couple of examples. 

First of all we describe the Painleve test applied to the equation ([3]), in 
which case y is just a scalar. Applying step 1 we look for leading order be- 
haviour which produces a singularity in the ODE, so we require y ^ cq (z—zq)'^ 
and /i must be a negative integer for a movable pole with no branching. This 
gives immediately 
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as the only possible dominant balance. Note that we could have also obtained 
this balance by assuming that y blows up as z zq, and then (since 52 is 
constant) >> 52 on the right hand side of the ODE, so the y^ term must 
balance with the left hand side of giving 

y"~62/^ as z ^ z^. (10) 

We can multiply by y' on both sides of pO|l and integrate to find 

i(2;')'^2y3, as z ^ zq (11) 

(throwing away the integration constant, which is strictly dominated by the 
other terms), and after taking a square root in (|11[) and integrating we find 
®. 

We now seek a solution of ([3]) given locally by a Laurent expansion around 
a double pole at z = zo, in the form 

00 

y = Y,Cj{z-zoy-\ co = l, (12) 
i=o 

where the value of cq has been fixed as in ([9]). We wish to determine the 
resonances^ which are the positions in the Laurent series p2p where arbitrary 
coefficients Cj can appear. Since the ODE ([3]) is second order, there must be 
two arbitrary constants in a local representation of the general solution: zq, 
the arbitrary position of the movable pole, and one other. To apply step 2 
of the Painleve test we take a perturbation of the leading order with small 
parameter e, in the form 

y^{z-zo)-^l + e{z~zor). (13) 
To first order in e we have 

2/2 ^ (z-zo)-4(l + 2e(z-zo)'-), y" ^ (z-zo)-'^(6 + e(r-2)(r-3)(z-zo)'^). 

Thus when we substitute the perturbation (|13p into the dominant terms (jlOp 
and retain only first order terms in e we find 



6y^ ^ e((r - 2)(r - 3) - I2) (z - z^y-^ = 0. 



Since the perturbation e is arbitrary, corresponding to the first appearance 
of a new arbitrary constant in the Laurent expansion (jl2p . the expression in 
large brackets must vanish, giving the resonance polynomial 



— 5r ~ 6 — 0, whence r = — 1 or r = 6. 
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The first resonance at r = — 1 must always be present in any expansion around 
a movable singularity, since it corresponds to the arbitrariness of zq. The sec- 
ond resonance value r = 6 indicates that the coefficient cg should be arbitrary. 

In order to complete the Painleve test, we must now substitute in the full 
Laurent expansion and check that it is consistent up to the coefficient cg. In 
this case we find that the expansion is precisely 

y = + ^3'^' - + 2^^^^" " + ■ ■ ■ ' 

so that cg = (73/28 is the arbitrary constant that appears after integrating ([3]) 
to obtain (j4]). In fact only even powers of (z — zq) occur in this expansion, since 
the Weierstrass function ([5]) is an even function of its argument. The higher 
coefficients in ()14|) can be found recursively in terms of the invariants 52, 33- 
(Up to overall multiples these coefficients are the Eisenstein series associated 
to the corresponding elliptic curve [91].) The pole position zq does not appear 
in the coefficients because the ODE ([3]) is autonomous. 

Here we should point out that passing the basic Painleve test is only a 
necessary condition for an ODE to have the Painleve property. Proving the 
Painleve property requires showing that the local Laurent expansions can be 
analytically continued globally to a single-valued function (or one with only 
fixed branched points), in the absence of movable essential singularities. For 
the ODE ([3]) this follows from the fact that the general solution ([5]) is given 
by a Weierstrass elliptic function, which is meromorphic (for a proof see e.g. 
[9TI I102j ). Painleve's proof that the first Painleve transcendent ([6]) is free 
from movable essential singularities is outlined by Ince in chapter 14 of [55] . 
but the proof is unclear and this has prompted recent efforts to find a more 
straightforward approach [HI jBH [H^ • 

Having seen an example where the Painleve test is passed, we now move on 
to an example for which it fails, by considering the following coupled second 
order system: 

y'{ = 2yl - 122/2, y'^ = 2yi2/2. (15) 

In |35j this system is associated to an interaction of four particles moving in a 
plane, subject to velocity-dependent forces, and in that context it is essential 
that both 2/1 (z), 2/2 (z;) (denoted C2(r), C4(r) in the original reference) and the 
independent variable z should be complex. To find the dominant balances, we 
look for leading order singular behaviour of the form 

2/1 ^ aZ^ 2/2 ~ (16) 

corresponding to a singularity in the solution at Z = z — Zq = for at least 
one of /i, V negative. Because the system ()15p is autonomous, we can expand 
in the variable Z, since the position zq of the movable singularity will not 
appear in the coefficients of local expansions around z = zq. 

There are three possible dominant balances for the system (fTSj) . namely 
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(i) yi ~ 3Z ^, 1/2 ^ bZ ^, b arbitrary; 
(a) Di ~ 3Z~^, ?;2 ^ bZ^ , b arbitrary; 

(m) yi ~ 10^-2, y2 f 

Other possible power law behaviour around Z = corresponds to /it, v both 
non-negative integers and leads to Taylor series expansions, which are not 
relevant to our analysis of singular points. 

The second step in applying the Painleve test is to find the resonances. For 
the system p3|) to possess the Painleve property we require that all resonances 
for all dominant balances must be integers, and at least one balance must have 
one resonance value of —1 with the rest being non-negative integers, in which 
case this is a principal balance for which the Laurent expansion should provide 
a local representation of the general solution. To find the resonance numbers 
r we substitute 

y^^aZ^'{l + 5Z''), y2 ^ bZ''{l + eZ"^) 

into the dominant terms of the system (jlSp for each of the balances (i) — (Hi), 
and take only the terms linear in S and e. This yields a pair of homogeneous 
linear equations for S, e (which correspond to the arbitrary coefficients ap- 
pearing at the resonances) . The determinant of this 2x2 system must vanish, 
which gives in each case a fourth order polynomial in r. 

Principal balance (i): It turns out that the balance (i) is the only principal 
balance, with resonances 

(i) r = -1,0,5,6. 

As mentioned before, the resonance —1 is always present, since it corresponds 
to the arbitrary position zq of the pole, while r = comes from the arbitrary 
constant b in the leading order term of the expansion for 1/2] the other two 
values arise from arbitrary coefficients higher up in the series for 7/1, 1/2, so 
that altogether there should be four arbitrary constants appearing in these 
Laurent series. However, for step 3 of the test we also require that all resonance 
conditions hold: so far we have only found the orders in the series where 
arbitrary constants may appear, but it is necessary to check that all other 
terms vanish at this order when the series are substituted into the equations. 
Taking 

00 00 
2/1 ^ L,{Z) := J2 KjZ', y2 - L2{Z) := ^ k2,jZ^ (17) 

i=-2 j=-2 

in the each of the equations (fT5|) we know already that the leading order terms 
require 

fci__2 — 3, ^2,-2 — b (arbitrary), 

giving the resonant term at r = in the expansion for y2 , while at subsequent 
orders we find 
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A:i _i = = fc2 -i; fci.o = ^2,0 = fci,i = = fc2,i- 

At the next orders we further obtain 

fci,2 = —36^/5, fc2,2 = 76'^/15; fci_3 = 0, /c2,3 arbitrary, 

so that the resonance condition at r = 5 corresponding to fc2,3 is satisfied. 
However, at the next order in the first equation of the system (If 5p . at the first 
appearance of the resonance coefficient fci,4, we find the additional relation 

fc2,2 = -6^5, 

which means that the resonance condition is not satisfied unless 6 = 0, con- 
tradicting the fact that b should be arbitrary. Thus the Painleve test is failed 
by this principal balance. 

The only way to rectify the failure of the resonance condition and leave b as 
a free parameter is to modify (fT7)l by adding logarithm terms. More precisely 
taking 

yi^Li(Z) + Zii(Z), L2iZ) + A^iZ), (f8) 

the resonance condition is resolved by taking 

Ai^-^b'zHogZ, A^r^-^b^ZHogZ. (f9) 

However, the additional terms Z\i, A2 in (fT8|) must then consist of a doubly 
infinite series in powers of Z and log Z , with the leading order behaviours 
given by (|f 9p. Only in this way is it possible to represent the general solution 
of the system (|f 5|) as an expansion in the neighbourhood of a singular point 
containing four arbitrary parameters. Such infinite logarithmic branching is a 
strong indicator of non- integrability [73l [87] . 

Non-principal balance (ii): The second balance denoted (m) above has 
resonances 

r = -5,-f,0,6. 

The presence of the negative integer value r = — 5 means that this is a non- 
principal balance. (For an extensive discussion of negative resonances see [20].) 
This gives Laurent expansions 

j/i - + fcZ-* - + O(Z^), y2^bZ^ + 0{Z'^). (20) 

In this case all resonance conditions are satisfied and all higher coefficients 
in (|20p are determined uniquely in terms of k and b. However, because it 
only contains three arbitrary constants (namely 5, k and the position zq of 
the pole), it cannot represent the general solution, but can correspond to a 
particular solution which is meromorphic. 

Non-principal balance (iii): For the balance {Hi) the resonances are given 
by r = — f and the roots of the cubic equation 
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- 15r2 + 26r + 280 = 0, 

which turn out to be a real irrational number and a complex conjugate pair, 
approximately 

r = -3.2676, 9.1338 ± 1.5048i. 

While non-integer rational resonances are allowed within the weak extension 
of the Painleve test (see [55] and section 6), irrational or complex resonances 
lead to infinite branching, and (as already evidenced by the principal balance 
(z)) the system ()15p cannot possess the Painleve property. This non-principal 
balance may be interpreted as a particular solution corresponding to a degen- 
erate limit of the general solution, and perturbation of this particular solution 
(within the framework of the Conte-Fordy-Pickering perturbative Painleve 
test [20]) will pick up the logarithmic branching present in the general solu- 
tion. Clearly it would have been sufficient to stop the test after the failure of 
the resonance condition in the principal balance (i), but we wanted to present 
the details of the other balances to show the different possibilities that can 
arise in the singularity analysis of ODEs. 

3 The Ablowitz-Ramani-Segur conjecture 

Having considered how to test for the Painleve property in ODEs, we now turn 
to the connection with integrable PDEs. In the 1970s it was discovered that 
ODEs of Painleve type, and in particular some of the Painleve transcendents, 
appeared as symmetry reductions of PDEs solvable by the inverse scattering 
technique. This led Ablowitz, Ramani and Segur [J to formulate the following. 

Ablowitz-Ramani-Segur conjecture: Every exact reduction of a PDE 
which is integrable (in the sense of being solvable by the inverse scattering 
transform ) yields an ODE with the Painleve property, possibly after a change 
of variables. 

To obtain ODE reductions of PDEs one can use the classical Lie symmetry 
method or its non-classical variants (see |77j for details) , or the direct method 
of Clarkson and Kruskal [16]. The idea is that having found the symmetry 
reductions of the PDE, one can either solve the ODEs that are obtained, or 
apply the Painleve test to them, to see if branching occurs. If all the ODE 
reductions are of Painleve type, then this suggests that the original PDE may 
be integrable. However, the need to allow for a possible change of variables 
will become apparent in section 6. Indeed, the most difficult aspect of this 
conjecture, if one would like to provide a proof of it, is in defining exactly 
what class of variable transformations should be allowed. 

As an example, consider the Korteweg~de Vries (KdV) equation for long 
waves on shallow water, which we write in the form 



Ut = Uxxx + QuUx- 



(21) 
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This has three essentially different reductions to ODEs; details of their deriva- 
tion are given in chapter 3 of |77| . The first is the travelling wave solution 

u{x,t)—w{z), z = X — ct, (22) 

where c is the (arbitrary) wave speed and w{z) satisfies 

w'" + 6ww' + cw' = 0. (23) 

After an integration and a shift in w this is equivalent to ([3|) , and the solution 
of (|23l) is given by 

w = -2p{z~ zo)-c/6, (24) 

where zg and the invariants 52 and 173 of the p-function are arbitrary constants. 
In the special case (72 = 4fc''/3, 53 = — 8fc^/27 the elliptic function degenerates 
to a hyperbolic function, and for c = — 4fc^ the reduction ([22)) yields the 
one-soliton solution 

u{x, t) = 2fc2sech^(fcx + Ak^t). (25) 

(Of course there is the additional freedom to shift the position of the soliton 
(|25p by the transformation a; — s- a; — xq.) 

The second reduction of KdV is the Galilean-invariant solution 

u{x,t) = -2{w{z) +t) , z^x-6f (26) 

where w{z) satisfies 

w'" - Uww' -1 = 0. (27) 

Upon integration, and making a shift in z to remove the constant of integra- 
tion, the ODE ([27]) becomes the first Painleve equation 

The third reduction of (PT|) is the scaling similarity solution 

u{x,t) = {-'Styiw{z), z = {-3ty^x. (28) 

This solution arises from the invariance of the PDE (|2ip under the group of 
scaling symmetries 

{x,t,u) — > (A.T, A^i, A^^u). 

After substituting the similarity form (|28|) into KdV and integrating once we 
find the ODE for w: 

w" + - zw\- = 0. (29) 

2w — z 

The parameter is the constant of integration, and ([29]) turns out to be 
equivalent to the equation P34, so called because it is labelled XXXIV in 
the Painleve classification of second order ODEs as detailed by Ince [SS] . The 
equation P34 can be solved in terms of the second Painleve equation ([7]), 
according to the relation 
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w = -y'-y'^, with £^a + l/2. (30) 

The above formula defines a Backlund transformation between the two equa- 
tions ([29]) and ([7]), and in fact there is a one-one correspondence between their 
solutions. With the parameters of the two ODEs related as in (|5n|) . the inverse 
of this transformation (defined for w ^ z/2) is given by 

w' + a 

y = o • 

For more details, and higher order analogues, see [49] and references. 

Thus we have seen that the ODE reductions of the KdV equation ([^ 
are solved either in elliptic functions or in terms of Painleve transcendents, 
and hence these reductions certainly have the Painleve property. So the KdV 
equation clearly fulfills the necessary condition for integrability required by 
the Ablowitz-Ramani-Segur conjecture, as it should do because it is inte- 
grable by means of the inverse scattering transform. In contrast to KdV, we 
consider another equation that models long waves in shallow water, namely 
the Benjamin-Bona-Mahoney (often referred to as BBM) equation fS], which 
takes the form 

Ut +Ux + UUx - Uxxt = 0. (31) 

The Bcnjamin-Bona-Mahoney equation is also known as the regularized long- 
wave equation, and was apparently first proposed by Peregrine [79J. The trav- 
elling wave reduction of the Benjamin-Bona-Mahoney equation is very similar 
to that for KdV: the PDE ^ has the solution 

u{x, t) = — 12cp(z — zq) + c — 1, z ^ X — ct, (32) 

given in terms of the Weierstrass p-function (with arbitrary values of the 
invariants g2, gs and the constant zq). In the hyperbolic limit with c = (1 — 
4fc^)~^ for k 7^ ±1/2 this gives the solitary wave solution 

X2fc2 

u(x, t) = ^ sech^(fca; - fc(l - 4fc2)"4), 

1 — 4fc^ 

but in contrast to ([25| this is not a soliton because the Benjamin-Bona- 
Mahoney equation is not integrable and collisions between such waves are 
inelastic: see the discussion and references in chapter 10 of [26] . 

Evidence for the non-integrable nature of the Benjamin-Bona-Mahoney 
equation is provided by another symmetry reduction, namely 

u{x^t) — -w{z) ~ \, z = X + nlogt^ (33) 

where k is a constant. Upon substitution of ([55]) into ([5T]) . w is found to satisfy 
the ODE 

Kw'" - w" - ww' - KW' + W = 0. (34) 
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For all values of the parameter k, the equation ([Mjl does not have the Painleve 
property, which means that (at least in these variables) the Benjamin-Bona- 
Mahoney equation fails the necessary condition required by the Ablowitz- 
Ramani-Segur conjecture. In the case k = 0, (j34p just becomes second order, 
so it is possible to compare with the list in chapter 14 of Ince's book [55] to 
see that w" + ww' — w = is not an ODE of Painleve type. A direct method, 
which works for any k, is to apply Painleve analysis directly to the equation 
and show that a resonance condition is failed. In fact the analysis can be 
greatly simplified by integrating in to obtain 

Kw" — w' KW = w(s) ds. (35) 

2 Jzi 

(The lower endpoint of integration zi is an arbitrary constant.) 

Now we can perform a Painleve test on the integro-differential equation 
(|35|) . For K ^ 0, in the neighbourhood of a movable singularity at z ^ zq the 
dominant balance is between the w" and w'^ terms, giving 

w{z) ~ V2k{z — 2;o)~^, z zq. 

If we suppose that this is the leading order in a Laurent expansion around 
z = zq, i-e. 

oo 

w{z) —'^^Wj{z — zqY''^ , = 12k, (36) 

i=0 

then at the next order we see that the coefficient of [z — ^o)^^ is 

Wl = - -, K 7^1/6. 

DK — 1 

(For K ~ l/Q the Laurent expansion immediately breaks down.) However, 
substituting the expansion ([36| into the left hand side of ([35|l gives a Laurent 
series, while on the right hand side there is a term log(z — zq) arising from the 
nonzero residue wi ^ Q. Hence the expansion (|36p cannot satisfy the equation 
([35]) . or equivalently (jM]) . and the Painleve test is failed. 

Thus we have seen that all of the ODE reductions of the KdV equation 
possess the Painleve property, but not all the reductions of the non-integrable 
Benjamin-Bona-Mahoney equation (pij) are of Painleve type. We leave it as an 
exercise for the reader to check whether the Benjamin-Bona-Mahoney equa- 
tion has other reductions apart from and ([551) (f^i' hints see exercise 3.2 
in [77]). However, it should be clear from the above that a fair amount of work 
is required when analysing a PDF in the light of the Ablowitz-Ramani-Segur 
conjecture, since one must first find all possible reductions to ODEs and then 
perform Painleve analysis on each of them separately. Finding the symmetry 
reductions can be a difficult enterprise in itself (see [17| for example), but in 
the next section we shall see how this complication can be avoided by using 
the direct method due to Weiss, Tabor and Carnevale [lOOj . 
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4 The Weiss- Tabor-Carnevale Painleve test 

While the symmetry reductions of PDEs are clearly indicative of their integra- 
bility or otherwise, it is more convenient to analyse the singularity structure of 
PDEs directly. This approach was pioneered by Weiss, Tabor and Carnevale 
|100| (hence it is usually referred to as the WTC Painleve test). However, in 
the context of PDEs with d independent (complex) variables zi, . . . ,Zd the 
singularities of the solution no longer occur at isolated points but rather on 
an analytic hypersurface S of codimension one, defined by an equation 



where (j) is analytic in the neighbourhood of S. The hypersurface where the 
singularities lie is known as the singular manifold, and it can be used to define 
a natural extension of the Painleve property for PDEs, which we state here 
in the form given by Ward [98] : 

Definition 4.1. The Painleve property for PDEs: If S is an analytic 
non- characteristic complex hypersurface in C^, then every solution of the PDE 
which is analytic on C'^\5 is meromorphic on C*. 

With the above definition in mind, it is natural to look for the solutions 
of the PDE in the form of a Laurent-type expansion near (/)(z) = 0: 



If the PDE has the Painleve property, then the leading order exponent 
appearing in the denominator of p8p should be a positive integer, with the 
expansion coefficients aj being analytic near the singular manifold = 0, and 
sufRciently many of these must be arbitrary functions together with the arbi- 
trary non-characteristic function (j). As mentioned in |60j in the context of the 
self-dual Yang-Mills equations, and further explained in j98j . it is important to 
state that </> should be non-characteristic because (even for linear equations) 
the solutions of PDEs can have arbitrary singularities along characteristics. 
The application of the Weiss- Tabor-Carnevale test using scries of the form 
proceeds as for the usual Painleve test for ODEs: when the series is sub- 
stituted into the PDE, equations arise at each order in (j) which determine 
the coefficients aj succesively, except at resonant values j = r, where the 
corresponding are required to be arbitrary (subject to compatibility con- 
ditions being satisfied). The Weiss- Tabor-Carnevale test is only passed if all 
resonance conditions are fulfilled for every possible balance in the PDE (i.e. all 
consistent choices of /i). Note that, just as for ODEs, passing the test merely 
constitutes a necessary condition for the Painleve property: a complete proof 
is much harder in general, although in the particular case of the self-dual 
Yang-Mills equations Ward [98] was able to use twistor methods to prove that 
they satisfy the requirements of Definition 4.1. 



<^(z)=0, {zi,...,Zd) eC^, 



(37) 




(38) 
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To see how the Weiss- Tabor-Carnevale test works, we wiU indicate the first 
steps of the analysis for the example of the KdV equation (PT|) . In that case, 
there are just two independent variables x and t, so d = 2, and there is only 
one dominant balance where the degree of the singularity for the linear term 
Uxxx matches that for the nonlinear term uux- Substituting an expansion of 
the form ([38|) into pT|) . with z = (x, i) € C^, it is clear that this gives /i = 2 
as the only possibility, and for the leading order and next-to-leading order the 
first two coefficients are determined as 

ao = -201, ai = 20J.J.. (39) 

This means that the expansion around the singular manifold for KdV can be 
written concisely as 

oo 

u{x,t) = 2{\og(j))xx + ^ak+2{x,t)(l){x,tt, (40) 

fc=0 

where it is necessary to assume so that (j) is non-characteristic. 

In general, at each order j there is a determining equation for the coeffi- 
cients of the series given by 

(j + l)(i - 4)(j - 6)a, = F.j[(j).,,(j)t,(j)xt ...,ak;k<j], (41) 

where the functions Fj depend only on the previous coefficients ak for k < j 
and their derivatives, as well as the various x and t derivatives of (j). It is clear 
from ()41l) that the resonance values are r = —1, 4, 6, meaning that we require 
0, a4 and ag to be arbitrary functions of x and t. For the KdV equation, 
apart from the standard resonance at —1 corresponding to the arbitrariness 
of (j), the other necessary conditions for r = 4, 6, namely = 0, Fg = 
are satisfied identically, and so in accordance with the Cauchy-Kowalevski 
theorem these three arbitrary functions are the correct number to provide 
a local representation (|40p for the general solution of the third order PDE 
([2T|) . We leave it to the reader to calculate the expressions for the higher Fj 
in (|4ip and verify the resonance conditions for F4 and F^; this is a standard 
calculation, so wc omit further details which can be found in several sources, 
e.g. (73i ilOO . For completeness we note that the issue of convergence of the 
expansion ([^D|) for KdV has also been completely resolved . 

We shall return briefly to the KdV equation in the next section, where 
we discuss how series such as (j38p can be truncated within the singular man- 
ifold method, leading to Backlund transformations and Lax pairs for inte- 
grable PDEs, and by further truncation to Hirota bilinear equations for the 
associated tau- functions. Before doing so, we would like to illustrate ways in 
which the basic Weiss- Tabor-Carnevale test may be further simplified, taking 
the non-integrable Benjamin-Bona-Mahoney equation (|3ip as our example. 
Applying the test as outlined above directly to the equation (pij) leads to 
an expansion pS]) very similar to that for KdV: it also has a single domi- 
nant balance with fi — 2 for a non-characteristic singular manifold (where 
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(jjx ^ ^ (j)t), and the same resonances r — —1,4,6, but for the Benjamin- 
Bona-Mahoney equation not aU resonance conditions are satisfied and the test 
is failed. It is a good exercise to perform this calculation and compare it with 
the corresponding results for KdV. Rather than presenting such a comparison 
here, we wish to give two shortcuts to the conclusion that the equation ((3T|) 
does not possess the Painleve property for PDEs. First of all, observe that if 
(j)x ^ then locally we can apply the implicit function theorem and solve the 
equation ([37| for x. Thus we set 

</)(x,i) =a;-/(i) (42) 

with f{t) df /dt ^ 0, and then we can take the coefficients in the expansion 
([55|) to be functions of t only; this is referred to as the 'reduced ansatz' of 
Kruskal, first suggested in [60]. With this ansatz, the Weiss- Tabor-Carnevale 
analysis for PDEs becomes only slightly more involved than applying the 
Painleve test for ODEs, and so constitutes a very effective way to decide if a 
PDE is likely to be integrable. 

For the Benjamin-Bona-Mahoney equation there is a second shortcut that 
can be made, which is to take the potential form of the equation by making use 
of the fact that it has a conservation law. This approach is widely applicable, 
since nearly all physically meaningful PDEs admit one or more conservation 
laws. For the equation (I3ip it is immediately apparent that it can be put in 
conservation form as 



du d f 1 , 



C = udx 



which implies that 



is a conserved quantity for the Benjamin-Bona-Mahoney equation, i.e. dC/dt — 
for u(x, t) defined on the whole real x-axis with vanishing boundary condi- 
tions at X = ±oo. It follows that upon introducing the potential v as the new 
dependent variable, with 



J — c 



u dx — > C as X ^ oo. 



we can replace u by w and its derivatives in (|3ip to obtain the potential form 
of the PDE, namely 

1 9 

vt - v.xxt + Vx + -vi^Q (43) 

(where we have integrated once and applied the boundary conditions to elim- 
inate the arbitrary function of t). If we now apply the Weiss- Tabor-Carnevale 
test to at the same time using the 'reduced ansatz' (UU, then we see 

that the only possible leading exponent in a Laurent-type expansion for v is 
/Li = 1, giving 
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vix,t)=Y./^At)i^-fit)Y''- (44) 

The equations for the coefficients Pj (t) at each order take the form 

(j + l)(j - - 6)/3, = F, [/, /,..., /3fe; A; < j], 

so the resonances are r = —1, 1,6 which compares with r = —1,4,6 for the 
original equation PT|) : clearly one of the resonances has shifted to a lower 
value by taking the equation in potential form (|43p . Upon substituting the 
series into the potential Benjamin-Bona-Mahoney equation, the leading 
order term is at order (f>~'^, giving the equation 

-6/3o/ + Ip^o = 0- 
Since /3o ^ 0, this determines the first coefficient as 

Po = 12/. 

However, at the next order (j>~^ in the equation ([43]) . we have the resonance 
r = 1 with the condition 

- 2/3o = 0, whence / = 0. (45) 

Since / is supposed to be an arbitrary non-constant function of t, we see that 
the resonance condition (|45|l is not satisfied, so the equation (|43|) fails the 
Weiss- Tabor-Carnevale Painleve test, indicating the non-integrability of the 
Benjamin-Bona-Mahoney equation. However, observe what happens if / is a 
linear function of t: then (|45)) is satisfied, corresponding to the travelling wave 
reduction ([32|l . which does have the Painleve property. 

The only way to remove the restriction (|45p on the function / would be to 
add a term —(/?//) log(x — f{t)) to the expansion It has been observed 
[55] that the inclusion of terms linear in log0 for PDEs in potential form is 
not incompatible with integrability. However, in this case terms of all powers 
of log(a; — f{t)) are required to ensure a consistent expansion in the potential 
Benjamin-Bona-Mahoney equation (I43p with three arbitrary functions /, /3i 
and Pq corresponding to the three resonances. 

For the reader who is interested in applying either the Painleve test for 
ODEs, as described in section 2, or the Weiss- Tabor-Carnevale Painleve test 
for PDEs, it is worth remarking that software implementations of these tests 
are now freely available. The web page www.mines.edu/f sJiome/whereman has 
algorithms written by D. Baldwin and W. Hereman, for instance. 

5 Truncation techniques 

Aside from the obvious application of the various Painleve tests in isolating 
potentially integrable equations (for example, in the classification of integrable 
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coupled KdV equations their usefulness can be extended by the means 

of truncation techniques. The first of these is known as the singular manifold 
method, which was primarily developed in a series of papers by Weiss |101j . 
The idea behind the method is that by truncating an expansion such as (I38p , 
usually at the zero order (0°) term, it is possible to obtain a Backlund trans- 
formation for the PDE. For such truncated expansions the singular manifold 
function (j) is no longer arbitrary, but satisfies constraints. In the case of in- 
tegrable equations that are solvable by the inverse scattering transform, the 
singular manifold method can be used to derive the associated Lax pair; for 
directly linearizable equations, such as Burger's equation or its hierarchy [80j . 
the method instead leads to the correct linearization. Even for non-integrable 
PDEs, where the constraints on cj) are much stronger, the singular manifold 
method can still be used to obtain exact solutions. Furthermore, for integrable 
PDEs the truncation approach can be carried further by cutting off the series 
before the zero order term, to yield tau-functions satisfying bilinear equations 

m- 

We will outline the basic truncation results for the KdV equation pTjl . 
before presenting more detailed calculations for the nonlinear Schrodinger 
(NLS) equation. For KdV, the Laurent-type expansion (j40|l can be consistently 
truncated at the zero order term to yield 

u = 2{\og(l))a;x+u, u = a2. (46) 

While substituting the full expansion ((40|l into KdV gives an infinite set of 
equations ((4T|l for (j) and the aj, the truncated expansion gives only a finite 
number. The last of these equations does not involve (j), and just says that u 
is also a solution of KdV, i.e. 

Ut = Uxxx + QuUx- 

The other equations (after some manipulation and integration) boil down to 
just two independent equations for and u, as follows: 



In the above, fc is a constant parameter. The important feature to note is 
that since u and u are both solutions of l[2T|) . the equation ((46|) constitutes a 
Backlund transformation for KdV, provided that cf) satisfies (|47p and For 
example, starting from the seed solution {t = 0, the Backlund transformation 
defined by ((46l) , ([iT]) and (|48l) can be used to generate the one-soliton solution 
((25|) . or even a mixed rational-solitonic solution by taking <j> — {x — \2kH) + 
{2k)-^smh{2kx + BkH). 



= 6fc^+(^-^)- (48) 
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It is maybe not immediately obvious that the system comprised of the two 
equations (|T7)) and (pS)) is equivalent to the standard Lax pair for KdV. This 
can be seen by making the squared eigenfunction substitution (l)^ — ijj^ , so 
that (j47p becomes a linear (time-independent) Schrodinger equation. In the 
context of quantum mechanics in one dimension, -0 is the wave function with 
potential —u and energy — fc^, i.e. (|47|) is equivalent to 

The second equation is known as the Schwarzian KdV equation [73] , and 
in its own right it constitutes a nonlinear integrable PDE for the dependent 
variable 4>; with the squared eigenfunction substitution it leads to the linear 
equation for the time evolution ipt- AH these results for KdV are well known, 
and have been extended to the whole KdV hierarchy; the interested reader 
who wishes to check these calculations is referred to [73] for more details. 

Perhaps less well understood, however, is the interesting connection [33] 
between the singularity structure of PDEs and the tau-function approach to 
soliton equations pioneered by Hirota [33 i which culminated in the Sato 
theory relating integrable systems to representations of afhne Lie algebras 
|7H I76j . The link with the singular manifold method is made by truncating 
the expansion (jlO]) at the last singular term in <j), and setting = r, to give 

w = 2(logT),„ (49) 

which is the standard substitution for the KdV variable u in terms of its tau- 
function. From (|^T|) . after substituting ((311) and performing an integration 
(subject to suitable boundary conditions), a bilinear equation is obtained for 
the new dependent variable r. This bilinear equation may be written concisely 
as 

(D^Dt - Dt)T • T = 0, (50) 
by making use of the Hirota derivatives: 

= ■ ^ - (i - ife) ' (s - 1^) ' 

The bilinear form is particularly convenient for calculating multi-soliton so- 
lutions ^46J , and leads to the connection with vertex operators JTTJ |7S1 [75] . 
For solitons the tau-function is just a polynomial in exponentials. In general 
r is holomorphic, so from it is clear that the places where r vanishes 
correspond to the singularities of u. 

We now present details on the application of the singular manifold method 
to the nonlinear Schrodinger equation 

«^t+V'xx-2|V'|V = 0. (51) 



This PDE (commonly referred to as NLS) describes the evolution of a complex 
wave amplitude 0, and due to the minus sign in front of the cubic nonlinear 
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term this is the non-focusing case of the nonhnear Schrodinger equation; the 
focusing case has +2|7/;p-0 instead, and describes a different physical con- 
text. The foUowing results on the singular manifold method for the nonlinear 
Schrodinger equation appeared in |47| . Seeking an expansion of the form (|38p 
for ([5T|) . at leading order we find the behaviour 

, aol =(t>x- 

9 

Thus, truncating the expansion at the zero order {cjf') level, we find 

V' = ^+^, ^ = ai- (52) 

To proceed with the singular manifold method we substitute the truncated 
expansion ([5^ into (|5ip. and set the terms at each order in (j> to zero. This 
yields the following four equations (the singular manifold equations): 

|aop-(/.2=0; 

: i(j)t + 2(t)^{\ogao)x + (j)^^ + 2ao?/' -I- 4ao'0 = 0; 
(j)^^ : iao,t + ao,xx - 4ao|^/'P - 2aoi/'^ = 0; 

0°: #t + ^,, -2|^|2V, = 0. 

Clearly the coefficient of (j>~^ just gives the leading order behaviour, while 
the cfP equation in ([53]) means that the truncated expansion ((52|) constitutes an 
auto-Backlund transformation for the nonlinear Schrodinger equation, since ip 
is another solution of ([5T|) . Observe that for x and t real, the singular manifold 
function (f) is seen to be real- valued from the leading order behaviour. Since the 
Painleve analysis is really concerned with singularities in the space of complex 
X, t variables, it is more consistent to write the nonlinear Schrodinger equation, 
together with its complex conjugate, as the system 

— — —2 (54) 
-iipt + ip^^-2^ ij; = 0, 

and then treat if) and tp as independent quantities. The system (|54p is the first 
non-trivial flow in the Ablowitz-Kaup-Newell-Segur (AKNS) hierarchy [3]. For 
this full system the singular manifold equations ([55)1 should be augmented 
with the corresponding 'conjugate' equations: formally these are obtained by 
taking the complex conjugate with 4> real (as for real x and t), and Oq, ip and 
ip complex. By formally taking the real and imaginary parts of the second 
equation in (|53p . which are equivalent to linear combinations of that equation 
together with its conjugate, the following consequences arise: 
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(t>xx + aQip + aoV^ = 0; ^^^^ 
i4't + (/)x(log[ao/ao])2: + SoV' - uqiI) = 0. 

Further manipulation of the singular manifold equations ()53p and their 
conjugates, together with ((55| . leads to the two equations 

ao,x = ~'2,iXao - 2tp4>^, (56) 

iao,t = (4A2 + 2?/'V;)ao + {-iiXi' + 2'iP^)<P^ (57) 

and their corresponding conjugates, where A is a constant. Upon substitution 
of the rearrangement 

ao = {ip- ■ip)(l) 



of ^ into ([56]), we find 

(V- - ^)x = -2zA(V' - V>) - (^ + V-) I - ^1 , (58) 

where we have used the first equation (j53p to substitute (/>x = |ao| = I^/j — V"! 
in the reduction to real x and t. A similar equation for (-0 — 4')t is obtained 
by eliminating ao and (p from (j57p , and the resulting relations between ip and 
tp together with (|58p constitute a Backlund transformation for the nonlinear 
Schrodinger equation in the form studied by Boiti and Pempinelli, taking the 
special case a = in the formulae of [10]. Starting from the vacuum solution 
tp = 0, and with zero Backlund parameter A = 0, this BT can be applied 
repeatedly to obtain a sequence of singular rational solutions of the nonlinear 
Schrodinger equation, which are described in [48] . 

The simplest singular rational solution has a single pole, which can be 
fixed at a; = 0. If we denote the sequence of these rational solutions {V'n}n>Oj 
then applying the BT ([55]) with A = starting from the vacuum solution the 
first three are 

with Ta being an arbitrary constant parameter which is real for real x and 
t. In general these rational functions can be written as a ratio of polynomial 
tau- functions 4'n = Gn/Fn satisfying bilinear equations (see below). The zeros 
and poles of each ■(/'„, which are the roots of the polynomials G„ and Fn re- 
spectively, evolve in t according to the equations of Calogero-Moser dynamical 
systems [48] . 

As well as leading to the Backlund transformation (|58p for the nonlinear 
Schrodinger equation, the singular manifold equations also yield the Lax pair, 
upon making the squared eigenfunction substitution 

"0 = -Xi, "0 = -xl- (60) 
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Fixing a sign we find immediately from tlie first equation (|53p that 

(f>x = XlX2, 

and then putting ([60]) into ([56|) . ([57]) and their conjugates gives a matrix 
system for the vector x = (xij Xa)^, that is 



Xx = Ux, 
Xt = Vx, 



(61) 



with the matrices 



(where we have replaced -0 by t/i in U, V). The system (j6ip is the non- focusing 
analogue of the Lax pair for the nonlinear Schrodinger equation found by 
Zakharov and Shabat [104| . and for U, V as above the PDE ([5T|) follows from 
the compatibility condition for the matrix system, which is the zero curvature 
equation 

Ut- V, + [U,V] =0. 

For real A, these matrices are elements of the Lie algebra su(l, 1), as opposed 
to su{2) for the case of the focusing nonhnear Schrodinger equation. 

To obtain the Hirota bilinear form of the nonlinear Schrodinger equation 
we can make a further truncation in (j52p . setting tp = Q, = G, (p — F ^ so 
that ((5T|) becomes 

^{{zDt + Dl)G ■ F) ~ ^{DlF ■ F + 2\G\^) = 0. 

The two equations in brackets can be consistently decoupled to give the bi- 
linear system for the two tau-functions F^ G: 



{iDt + Dl)G ■ F ^ 0; 
DlF-F + 2\G\^ ^0. 



(62) 



It is easy to check that the numerators and denominators in the rational 
functions ([5^ are particular solutions of the system ([5^ . The bilinear form 
of the nonlinear Schrodinger equation was used by Hirota to derive compact 
expressions for the multi-soliton solutions [JS]. A further consequence of (p^ 
is the bilinear equation 

iD^DtF ■ F -2D^G -G = i^F^, (63) 

with a constant 7. This constant can be removed by a gauge transformation 
of the tau-functions, rescaling both F and G by exp[7a;t/2]. Eliminating G 
between (|63)) and ([62]), the nonlinear Schrodinger equation is then rewritten 
as a single trilinear equation, expressed as a sum of two determinants, namely 
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The tau-function solution of the trihnear equation is sufBcicnt to deter- 
mine both the modulus and the argument of the complex amplitude ip (see 
[48] and references). 

From the preceding results for the KdV and nonlinear Schrodinger equa- 
tions it should be clear that truncation methods can be extremely powerful in 
extracting information about integrable PDEs. There are several refinements 
of the singular manifold method, in particular those involving truncations us- 
ing Mobius-invariant combinations of (j) and its derivatives [HI [72], and the use 
of two singular manifolds for PDEs with two different leading order behaviours 
[21] . Probably the most elegant and general synthesis of these extended meth- 
ods is the approach formulated by Pickering ^81j, who uses expansions in a 
modified variable satisfying a system of Riccati equations. Truncation methods 
have even been used to derive Backlund transformations for ODEs, in partic- 
ular Painleve equations [18] . However, it is uncertain whether such methods 
can really be made sufficiently general in order to constitute an algorithmic 
procedure for deriving Lax pairs for integrable sytems. In particular, trun- 
cation methods are not directly applicable to integrable PDEs which exhibit 
movable algebraic branching in their solutions, which are the subject of the 
next section. 



6 Weak Painleve tests 

There are numerous examples of integrable systems which do not have the 
strong Painleve property, but which satisfy the weaker criterion that their 
general solution has at worst movable algebraic branching. Perhaps the sim- 
plest example is to consider a Hamiltonian system with one degree of freedom 
defined by the Hamiltonian (total energy) 

H^\p' + V{q), 

where the potential energy V is a polyomial in q of degree d> 5. The equations 
of motion (Hamilton's equations) are 

dq dp 
which are trivially integrable by a quadrature: 
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If the potential energy is normalized so that the leading term of the polynomial 
is —2q'^/{d — 2)^, then with q{t) having a singularity at t = to the integral in 
(|55|) gives 

, n {2^d)dQ 1 

^^l,"^ =±q'-"^ as q^^. 

(for a suitable choice of branch in the square root). Thus at leading order we 
have 

q^±{t~t,f/(^-^'). (66) 

For both d = 2g + 1 (odd) and d = 2g + 2 (even) q is determined by the 
hyperelliptic integral (|65p corresponding to an algebraic curve of genus g. 
When (7=1 the solution is given in terms of Weierstrass or Jacobi elliptic 
functions, and both q and p are meromorphic functions of t. However, for a 
potential of degree 5 or more we have g > 2, and it is clear from that 
q has an algebraic branch point at t ~ Iq, since in that case 2/(2 — d) is a 
non-integer, negative rational number. In fact it is easy to verify that l|66p is 
the leading order term of an expansion in powers of {t - to)^/<-'^-^\ Rather 
than being meromorphic as in the elliptic case, for d > 5 the function q(t) 
is generically single-valued only on a covering of the complex t-plane with 
an infinite number of sheets, and has an infinite number of algebraic branch 
points (see [2]). 

Clearly for potentials of degree 5 or more, this simple Hamiltonian system 
fails the basic Painleve test, and yet it is certainly integrable according to any 
reasonable definition. (Indeed, any Hamiltonian system with one degree of 
freedom is integrable in the sense that Liouville's theorem holds.) In order to 
avoid excluding such basic integrable systems from singularity classification, 
Ramani et al. [86] proposed an extension of the Painleve property. 

Definition 6.1. The weak Painleve property: An ODE has the weak 
Painleve property if all movable singularities of the general solution have only 
a finite number of branches. 

There are many examples of finite-dimensional many-body Hamiltonian 
systems which are Liouville integrable and yet have algebraic branching in 
their solutions [H [2] • Among these examples [2] is the geodesic flow on an 
ellipsoid, which was solved classically by Jacobi [58]. Many other examples, 
such as those considered by Abenda and Fedorov in [1], arise naturally as 
stationary or travelling wave reductions of PDEs derived from Lax pairs, in 
particular those obtained from energy-dependent Schrodinger operators '50|. 
Thus the corresponding Lax-integrable PDEs have algebraic branching in their 
solutions, and fail the Weiss- Tabor-Carnevale test described in section 4. It is 
natural to extend the notion of the weak Painleve property to PDEs as well, 
and perform Painleve analysis on ODEs and PDEs with this property by 
allowing algebraic branching and rational (not necessarily integer) values for 
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the resonances. We illustrate this procedure with the example of the Camassa- 
Holm equation and a related family of PDEs 51j which have peaked solitons 
(peakons). 

The Camassa-Holm equation was derived in |12j by asymptotic methods as 
an approximation to Euler's equation for shallow water waves, and was shown 
to be an integrable equation with an associated Lax pair. In the special case 
when the linear dispersion terms are removed the equation takes the form 

Ut - Uxxt + iuUx ^ 2UxUxx + UUxxx, (67) 

and in this dispersionless limit it admits a weak solution known as a peakon, 
which has the form 

u(a;,i) = ce-l^^-"*!. (68) 

Note that the notion of a 'weak solution' (as defined in [37], for instance) is 
completely unrelated to the 'weak' Painleve property. The peakon solution has 
a discontinuous derivative at the position of the peak, and the dispersionless 
Camassa-Holm equation (|67|) has exact solutions given by a superposition of 
an arbitrary number of peakons which interact and scatter elastically, just as 
for ordinary solitons. A detailed analysis of weak solutions of ([57)1 has been 
performed by Li and Olver i68j. 

However, the Camassa-Holm equation is an example of an integrable equa- 
tion which does not satisfy the requirements of Definition 4.1, but instead 
passes the weak Painleve test. In the neighbourhood of an arbitrary non- 
characteristic hypersurface (f>{x,t) — where the derivatives of u blow up, it 
admits an expansion with algebraic branching: 

oc 

u{x,t) = -(j>t/c^x + J2a,{x,t)c^^/^+^/\ (69) 
i=o 

If we regard the branching part (jy^^^^ as the leading term (since it produces 
the singularity in the derivatives m^, Ut on = 0), then the resonances are 
r = —1,0,2/3 which correspond to the functions (/>, ao, a2 being arbitrary. The 
Camassa-Holm equation thus satisfies the weak extension of the Weiss- Tabor- 
Carnevale test, since the expansion (|69p is consistent, with the resonance 
conditions at r = and r = 2/3 being satisfied. Of course the test is only 
local, whereas the weak Painleve property is a global phenomenon, and to 
prove it rigorously for this PDE would require considerable further analysis. 
The weak extension of the Painleve test is still a useful tool, in the sense that 
if an equation has irrational or complex branching (either at leading order 
or in its resonances), or if a failed resonance condition introduces logarithmic 
branching into the general solution, then this is a good indication of non- 
integrability. Nevertheless, even for ODEs the weak Painleve property should 
be applied cautiously as an integrability criterion. For an excellent discussion 
see HI]. 
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We would now like to apply the weak Painleve test to a one-parameter 
family of PDEs that includes (|67p . before showing the effect that changes of 
variables can have on singularity structure. We shall consider the family of 
PDEs 

Ut - Uxxt + {b+ l)uUx ^ hUxUxx + UUxxx, (70) 

where the parameter b is constant. These are non-evolutionary PDEs: due 
to the presence of the Uxxt term, (|70p is not an evolution equation for u. 
The (dispersionless) Camassa-Holm equation is the particular member of this 
family corresponding to 6 = 2. The original reason for interest in this family 
is that Degasperis and Procesi applied the method of asymptotic integrability 
[24j and isolated a new equation as satisfying the necessary conditions for 
integrability up to some order in a multiple-scales expansion. After removing 
the dispersion terms by combining a Galilean transformation with a shift in 
u and rescaling, the Degaspcris-Procesi equation can be written as 

Ut - Uxxt + '^UUx = iUxUxx + UUxxx, (71) 

which is the & = 3 case of ((TO)) , and it was proved in [25] by construction of 
the Lax pair that this new equation is integrable. A powerful perturbative 
extension of the symmetry approach was also applied to the non-evolutionary 
PDEs ([70]) in [70], and it was confirmed that only the special cases h — 2 
(Camassa-Holm) and 6 = 3 (Degasperis-Procesi) fulfill the necessary condi- 
tions to be integrable. Hamiltonian structures and the Wahlquist-Estabrook 
prolongation algebra method for these PDEs have also been treated in detail 
[51j . Subsequently it has been shown that (after including dispersion) every 
member of the family ((70]) arises as a shallow water wave equation [25 , except 
for the special case 6 = — 1. 

For Painleve analysis it is convenient to rewrite (|70p in the form 

rrit + urux + buxTTi — 0, m = u — Uxx- (72) 

To apply the weak Painleve test, we look for algebraic branching similar to 
the leading order in ()69|) . with the derivatives of u blowing up on a singular 
manifold <t>{x,t) — 0. Thus we seek the following leading behaviour: 

u^UQ + a<j)^, /leZ, 0</i<l. (73) 

Then for the derivatives of u and m as defined in (|72p the most singular terms 
are as follows: 

Ux ^ a4>xtJ.(t>^~^ ■, m -a4>lfi{p - 1)4>^~^, 

rux - -c^4>lf^{p - 1)(m - 2)(/)""^ mt - -afytKfJ- " 1)(m " 2)0''-^ 

Substituting these leading orders into (|72p we find a balance at order 0^"^ 
between the nit and unix terms provided that 
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"0 = -<t>tl<t>x- 

The next most singular term in the PDE is then at order 0^'^"^, corresponding 
to a balance between the unix and Uxm terms in (1721) . with coefficient 



and this is required to vanish giving 

Thus we see that for a weak Painleve expansion with the leading exponent 
H being a rational number between zero and one, the most singular terms 
require that the parameter h should also be rational with 

2 

6 = - - 1 > 1. 

To find and test the resonances in an expansion with this leading order, it is 
sufficient to take the reduced ansatz p2)) for and then make a perturbation 
of the leading order terms with parameter e: 

j{t)+a{tW{l + e<f ), cf>^x~f{t). (75) 



Substituting the perturbed expression into (|72[) and keeping only terms linear 
in e, we see that terms possibly appearing at order (/)^+''^'^ cancel out auto- 
matically (due to the form of uq), leaving the resonance equation coming from 
the coefficient of which is 

-ea^{r^ + (2/i - l)r^ + 2(/i - l)r) = 0. 

Hence the resonances are 



r = -l, 0, 2(1 -Ai), 

with given in terms of the parameter b by (|74p . 

Having applied the first part of the weak Painleve test and found a dom- 
inant balance and the corresponding values for the resonances, it becomes 
apparent that the test is completely ineffective as a means to isolate the two 
integrable cases b — 2 and 6 = 3 of ([7^ . Although the leading order reso- 
nance r = (corresponding to a being arbitrary) is automatically satisfied, 
the second resonance condition at r = 2(1 — /x) must be checked for every 
rational value of /i with < /i < 1 (or equivalently every rational value of 
the parameter b > 1). If we write in its lowest terms as a ratio of positive 
integers, /i — N1/N2, then l|73p is the leading part of an expansion for u in all 
powers of (/)^/^^ , and as the difference N2 — Ni increases there is an increas- 
ingly large number of terms to compute before the final resonance is reached. 
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Checking this resonance for the whole countable infinity of rational numbers 
6 > 1 seems to be a totally intractable task. Gilson and Pickering showed that 
all the PDEs within a class including ([7^ failed every one of a combination 
of strong Painleve tests [Mj. Nevertheless, it is simple to verify that the weak 
Painleve test is satisfied for the two particular cases 5 = 2,3 which are known 
to be integrable. 

However, after a judicious change of variables, involving a transformation 
of hodograph type, it is still possible to use Painleve analysis to isolate the 
two integrable peakon equations. Such transformations have been applied to 
integrable PDEs with algebraic branching (see [15]) in order to obtain equiva- 
lent systems with the strong Painleve property. That this should be possible is 
in accordance with the Ablowitz-Ramani-Segur conjecture, but the difficulty 
lies in finding the correct change of variables. In fact, for a general class of 
systems that display weak Painleve behaviour (related to energy-dependent 
Schrodinger operators) we presented a particular transformation in [50| and, 
from an examination of a principal balance, we asserted (without proof) that 
this transformation produced equivalent systems with the strong Painleve 
property. However, from a more careful calculation of other balances we have 
recently observed that this earlier assertion was incorrect [53]. In the case of 
the Camassa-Holm equation ([57]) . a link to the first negative flow in the KdV 
hierarchy was found by Fuchssteiner [32], and in [51] it was shown that the 
appropriate transformation can be extended to (almost) every member of the 
family of non-evolutionary PDEs (f72|) . 

The key to a suitable change of variables for (|72)) is the fact that for any 
is a conserved quantity, with the conservation law 

Pt -{pu)x, m = -p^. (76) 

This allows a reciprocal transformation, defining new independent variables 
X, T via 

dX — pdx — pu dt, dT ~ dt. (77) 

Observe that the closure condition d^X — Q for the exact one- form dX is 
precisely ()76|) . and transforming the derivatives yields the new conservation 
law 

{p'^)t=ux. (78) 
In the old variables, p is related to u by 

/ - {dl - l)u, (79) 

Replacing dx by pdx and using (l78]) . this means that ([79]) can be solved for 
u to give the identity 

-p(logp)xT -p''- (80) 

Finally the conservation law (|78|) can be written as an equation for p alone, 
by substituting back for u as in (I80|) to obtain 



Painleve tests, singularity structure and integrability 29 

Thus we have seen that for each 6 7^ 0, the equation ((72|) is reciprocally 
transformed to (|8ip , with the new dependent variable p and new independent 
variables X, T as in ([77)1 . (For more background on reciprocal transforma- 
tions, see [64].) By making the subsitution p — exp{iri), ((8T|) becomes a gener- 
alized equation of sine-Gordon type |51j . The point of making the reciprocal 
transformation is that we may now apply the strong Weiss- Tabor-Carnevale 
Painleve test to the equation in these new variables. At leading order near a 
hypersurface 4i{X, T) = there are two types of singularity that can occur in 
the equation (I8ip , corresponding to p either vanishing or blowing up there: 

• p ^ acj), for b> —1, with a = 'i^4>^ for 6 7^ —1; 

• p~ /30^, for = 2/(1 - 6) < 1. 

In the first balance, the resonances are r — —1, 1, 2. However, if we require the 
strong Painleve test to hold we see that we must have 6 S Z, since otherwise 
the p^ term will introduce branching into the expansion in powers of 0. The 
second balance can only hold for |6| > 1, but if < —1 then /i ^ Z, while if 
6 > 1 then requiring /i = 1 — M to be a (negative) integer gives 

^=^7^' M = 2,3,4,... (82) 

From the first balance we require b to be an integer, and the only integer values 
in the sequence ([5^ are 5 = 2,3 (corresponding to M = 3, 2 respectively). 
Interestingly, when the Wahlquist-Estabrook method is applied to ((7^ . this 
same sequence crops up from purely algebraic considerations |51| . 

The above analysis shows that the two integrable cases 6 = 2,3 are isolated 
immediately just by looking at the leading order behaviour. It is then straight- 
forward to show that for both types of singularity in the equation ([5T|) , these 
two cases fulfill the resonance conditions and thus satisfy the strong Painleve 
test. However, the observant reader will notice that further analysis is required 
to exclude the two special integer values 6 = ±1, for which only the first type 
of singularity arises; this is left as a challenge to the reader. 



7 Outlook 

It should be apparent from our dicussion that the various Painleve tests are 
excellent heuristic tools for identifying whether a given system of differen- 
tial equations is likely to be integrable or not. However, the strong Painleve 
property is clearly too stringent a requirement, since it is not satisfied by a 
large class of integrable systems which have movable algebraic branch points 
in their solutions. On the other hand, checking all possible resonances in the 



30 Andrew N.W. Hone 



weak Painleve test can be impractical as a means to isolate integrable sys- 
tems, and if there are negative resonances then more detailed analysis may 
be necessary to pick up logarithmic branching '82]. In this short review we 
have concentrated on methods for detecting movable poles and branch points. 
However, for equations like (|67|). the existence of the peakon solution ([68]) has 
led to the promising suggestion that Dirichlet series (sums of exponentials) 
may be a useful means of testing PDEs [84j- Also, although we have only con- 
sidered singularities of ODEs in the finite complex plane, there are extensive 
techniques for analysing asymptotic behaviour at infinity [501 IMl ISH] • 

Before closing, we should like to give a brief mention to the fruitful connec- 
tion between the singularity structure and integrability of discrete systems, 
in the context of birational maps or difference equations. In the last twenty 
years, there has been increased interest in discrete integrable systems. Liou- 
ville's theorem on integrable Hamiltonian systems extends naturally to the 
setting of symplectic maps or more generally to Poisson maps or correspon- 
dences 197] , and many new examples of integrable maps have been found 
[94] . Grammaticos, Ramani and Papageorgiou introduced a notion of singu- 
larity confinement for maps or difference equations I38j , which they used very 
successfully as a criterion to identify discrete analogues of the Painleve equa- 
tions, and they proposed that it should be regarded as a discrete version of 
the Painleve property. 

In order to illustrate singularity confinement, we shall consider the second 
order discrete equation 

u„+i(u„)^u„-i = aq^Un -f /?, (83) 

which is a non-autonomous version of an equation of the Quispcl-Roberts- 
Thompson type [85] . and can be explicitly solved in elliptic functions in the 
autonomous case q = I [52] ■ For g ^ 1 the equation (|83]) can be regarded 
as a discrete analogue of the first Painleve equation, because if we set Un = 
h"^ — y{nh), a — 4/i~®, (3 — — 3/i~®, q — 1 ~ /A and take the continuum 
limit /i — > 0, with z = nh held fixed, then equation ^ arises at leading order 
in h. 

The idea of singularity confinement is that if a singularity is reached upon 
iteration of a discrete equation or map, then it is possible to analytically 
continue through it. (This is by analogy with the fact that the solution of an 
ODE with the Painleve property has a unique analytic continuation around 
a movable pole.) In the case of (|83]) . a singularity will be reached if one of 
the iterates, say un, is zero, because this means that the next iterate un+i is 
not defined. By redefining a and shifting the index n if necessary, we can take 
A'' = 1 without loss of generality, so ui = 0. The vanishing of ui requires that 
at the previous stage auo + P = must hold. Setting u-i = a (arbitrary) and 

auQ + (3 — e 

gives Ml ~ P^^a^^e — > as e — > 0, and the singularity appears at 
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U2 ~ -fi^a'a-^e'^ 

However, subsequently we have ~ —q^a^ I3~'^a~^e, U4 = 0(1) and further 
iterates are regular in the limit e — s- 0. In this sense, we say that the singularity 
is confined. 

Although the singularity confinement criterion led to the discovery of many 
new discrete integrable systems (see j88j and references), it was shown by 
Hietarinta and Viallet that it is not a sufficient condition for integrability 
|42| . In fact, they found numerous examples of maps of the plane defined by 
difference equations of the form 

Un+l + Un-1 = fiUn), 

for certain rational functions /, which have confined singularities and yet 
whose orbit structure displays the characteristics of chaos. Other examples 
of singularity confinement in non-integrable maps can be found in [54] . Nev- 
ertheless, it seems that singularity confinement should be a necessary condi- 
tion for integrability of a suitably restricted class of maps. In fact, Lafortune 
and Goriely have shown that for birational maps in d dimensions, singularity 
confinement is a necessary condition for the existence of d — 1 independent 
first integrals [S7]. Ablowitz, Halburd and Herbst have made an alternative 
proposal for extending the Painleve property to difference equations, by using 
Nevanlinna theory [7l|4T], and this has deep connections with various algebraic 
or arithmetic measures of complexity in discrete dynamics (see [40l |42l \89\ l92] 
and references). 

For the reader who is interested in pursuing the subject of Painleve analysis 
and its applications to both integrable and non-integrable equations, a number 
of excellent review articles are to be recommended [291 ESI ESI [87l [95] , as well 
as the proceedings volume [22] . 
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